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ABSTRACT 

A numerical solution to the problem of hydrodynamic lubrication of rigid 
point contacts with an isoviscous, incompressible lubricant has been obtained. 
The hydrodynamic load-carrying capacity under unsteady (or dynamic) conditions 
arising from the combined effects of squeeze motion superposed upon the en- 
training motion has been determined for both normel approach and separation. 
Superposed normal motion considerably increases net load-carrying capacity 
during normal approach and substantially reduces net load-carrying capacity 
during separation. Geometry has also been found to have a significant influ- 
ence on the dynamic load-carrying capacity. The ratio of dynamic to steady 
state load-carrying capacity increases with increasing geometry parameter for 
normal approach and decreases during separation. The cavitation (film rupture) 
boundary is also influenced significantly by the normal motion, moving down- 
stream during approach and upstream during separation. For sufficiently high 
normal separation velocity the rupture boundary may even move upstream of the 
minimum-film-thickness position. 

*6anaras Hindu University, Varanasi, India, and NRC-NASA Research Associate. 
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sixty-three cases were used to derive a functional relationship for the 

ratio of the dynamic to steady state load-carrying capacity e in terms of 

the dimensionless normal velocity parameter q (incorporating normal velocity, 

entraining velocity, and film thickness) and the geometry parameter a. The 

result is expressed in the form 

0 ja” 0'028 gg qjj 1/q 

The ratio of the dynamic to steady state peak pressures in the contact 
5 increases considerably with increasing normal velocity parameter during 
normal approach, with a similar decrease during separation. The ratio is ex- 
pressed as a function of q and o by 

5 = sech (2q)| 
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NOMENCLATURE 

Amplitude of harmonic motion, m 

Film thickness, m 

Fluid inlet level . m 

Central (or minimum) film thickness, m 

Dimensionless film thickness, h/R^ 

Dimensionless fluid inlet level, h^_/R„ 

Dimensionless minimum (or central) film thickness, hg/R^ 

Pressure, N/m^ 

Maximum or peak pressure along line of minimum film thickness, 
N/m^ 

Dimensionless pressure, pRj^/ngUg 

Dimensionless peak pressure along line of minimum film thickness, 
N/m^ 

1/2 

Dimensionless normal velocity parameter, Uj^j/Ug ( 1/2 Hq) 

Effective radius of curvature, RyRJ(R^ 

X y X y 

Separation due to geometry of solid, m 

Surface velocity of solids A and B 

Normal velocity of approach or separation of solids, m 


( 




L 


Dimensionless normal velocity^ 

Uj Average surface velocity, (U^ + U0)/2, m 

t Time, sec 

X Coordinate along rolling direction, m 

Xp Location of peak pressure from minimum-film-thickness position, m 

x^ Location of film rupture boundary from minimum-film-thickness 

position, m 

X Dimensionless coordinate, x/R^^ 

Xp Dimensionless location of peak pressure from minimum-film- 

thickness position on line of minimum film thickness, Xp/R^^ 

Xj, Dimensionless location of film rupture boundary from minimum-film- 

thickness position on line of minimum film thickness, x^/R^^ 

y Coordinate transverse to rolling direction, m 

Y Dimensionless coordinate, ylR^ 

w Load-carrying capacity, N 

W Dimensionless load-carrying capacity, w/tiqU^Rj^ 

a Radius ratio (geometry parameter), Ry/R^^ 

a Dynamic load ratio (ratio of dimensionless dynamic to steady 

state load-carrying capacity), W/(W)q = 0 

e Amplitude ratio, a/(h«) 

^ mean 

2 

Dq Fluid viscosity at standard temperature and pressure, N s/m 

5 Dynamic peak pressure ratio (ratio of dimensionless dynamic to 

steady state peak pressure, P^jax^^^^max^ q = 0 
to Circular frequency of harmonic oscillation 

INTRODUCTION 

Hydrodynamic lubrication of rigid line and point contacts has attracted 
the attention of tribologists for the past two decades because of its practical 
applications in the lubrication of lightly loaded rolling-element bearings. 

The major focus has been in estimating the minimum film thickness and in ob- 
taining a relationship between film thickness, load, and entraining velocity 
parameters. An early theoretical approach to the hydrodynamic lubrication of 
rigid cylinders for an incompressible, isoviscous lubricant was presented by 
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Martin (1916). A classical analytical solution for rigid spherical contacts 
using half-Swnmerfeld boundary conditions for the film rupture was obtained by 
Kapitsa (1955). Experimental measurements of the film thickness in hydro- 
dynamically lubricated point contacts were reported by Dalmaz and Godet (1974) 
and Thorp and Gohar (1972). Brewe et al. (1979) obtained a film thickness 
equation for the lubrication of fully flooded, rigid, isoviscous point contacts 
through a numerical analysis that used a more realistic Reynolds boundary con- 
dition for the film rupture In the exit region. They showed that the film 
rupture boundary condition resulted in an additional geometry parameter in the 
film thickness equation. Inlet starvation effects have also been irwestigated 
by several researchers, viz Chiu (1974), Dowsor (1968), Floberg (1965, 1973), 
and Wedeven et a1. (1971). More recently inlet starvation effects were incor- 
porated into the film thickness equation by Brewe and Hamrock (1982). 

Sasaki et al. (1962) presented a solution for the isothermal lubrication 
of rigid cylinders subjected to sinusoidal load fcr a non-Newtonian fluid by 
using superposition of pressure curves generated by normal approach and en- 
training velocities. Appropriate boundary conditions for the film rupture in 
the exit were not satisfied. Vichard (1971) analyzed theoretically the tran- 
sient effect associated with squeeze film action under both hydrodynamic and 
elastohydrodynamic conditions. He concluded that at low film thickness the 
damping phenomenon associates' with normal approach was more important under 
elastohydrodynamic conditions than under hydrodynamic conditions. Investiga- 
tions of the elastohydrodynamic theory of concentrated contacts in normal 
approach were carried out by Herrebrugh (1970), Christensen (1962), and Lee 
and Cheng (1973), and more recently by Wang and Cheng (1981) and Chandra and 
Rogers (1983), 
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Dowson, Markho, and Jones (1976) presented a general theoretical analysis 
of the hydrodynamic lubrication of rigid cylindrical contacts by an isoviscous 
lubricant in combined rolling and normal motion. Results showed that normal 
motion significantly influences the load-carrying capacity of the contact. 

The film rupture boundary was also observed to be affected significantly by 
normal motion. Experimental investigations were also reported by Markho and 
Dowson (1976). However, a similar theoretical analysis does not exist that 
incorporates the geometry effect into the hydrodynamic lubrication of rigid 
point contacts by an isoviscous, incompressible lubricant in combined rolling 
and normal motion. The present investigation deals with this problem. Through 
a numerical analysis the effect of squeeze film action due to normal motion 
combined with entraining action due to pure rolling has been investigated for 
rigid point contacts. Geometry effects have also been incorporated. Computer 
plots of pressure distribution depict the effect of normal motion on the pres- 
sure distribution and film rupture boundary in the exit region. Sixty-three 
cases were used to derive an empirical formula for the dynamic load-carrying 
capacity as a function of the dimensionless normal velocity and geometry pa- 
rameters. The dynamic peak pressure ratio is also given by a simple formula 
as a function of these parameters. 


THEORETICAL ANALYSIS 


Reynolds Equation 

The Reynolds equation for the hydrodynamic lubrication of two rigid solids 
separated by an incompressible, isoviscous lubricant film is given as Eq. (1). 
Lubricant flow for a rolling-sliding contact and the corresponding pressure 
buildup is shown in Fig. 1. 



12r)oUs 
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where 

X * x/R^; Y « y/R^i H * h/R^; P « pRu/r^U^; a « R^/R^; 

Uj = 1/2 (U^ + Ug); = Uj^/Ug 

Eq. (P can then be expressed In dimensionless form as 


( 2 ) 



where q « Uj^/ is the dimensionless normal velocity parameter that 

incorporates the dimensionless normal velocity and the central film 
thickness Hq. 


The boundary conditions for Eq. (4) are given as 
P = 0 at the inlet boundary (i.e., H = \ 

dP I 

P = ^ = 0 at the cavitation boundary ) (5) 

(i.e., Reynolds boundary l 

condition), where N is the normal ) 
to the boundary. ' 


F i Im Shape 

The thickness of a hydrodynamic film between two rigid bodies in rolling 
contact can be written as the sum of two terms, that is, 

h = ho + S(x,y) (6) 

where 

hg Central (or minimum) film thickness 
S(x,y) Separation due to geometry of solids 
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The separation of two rigid solids is shown in Fig. 2. The film shape 
can be expressed in dimensionless form by following the widely used parabolic 
approximation as 

u2 y2 

+ (7) 

The dimensionless parameters involved in the subsequent numerical investigation 
are therefore the central film thickness Hq, the normal velocity parameter 
q, and the geometry parameter a. 

A pressure distribution satisfying the Reynolds equation (Eq. (4)) and 
the boundary conditions given by Eq. (5) was determined for a given speed, 
vis':osity, central film thickness, fluid inlet level, and normal velocity par- 
ameter by using the Gauss-Seidel Iterative method with overrelaxation. For 
optimum efficiency nonuniform nodal structure consisting of coarse and fine 
mesh was used to enhance accuracy in the region of high pressure and large 
pressure gradients. 

A fine mesh spacing of 0.001 and a coarse mesh spacing of 0.01 were used 

in all of the cases run for a =1.0. For geometry parameters other than 

a a 1.0 boundaries of the computation zone were located according to the rela- 
0 *5 

tionship Y = a X and correspondingly coarser grids in the Y direction 
were used for higher values of a. Grid sizes were varied until no appreciable 
influence on the pressure distribution or the load-carrying capacity was found. 

Hydrodyanmic Load-Carrying Capacity 

The load-carrying capacity can be calculated by integrating the pressure 
in the contact region and is written as 

w = J p dxdy (8) 

A 
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In dimensionless form it is expressed as 

W - w/r^UgR^ * / P dxdy (9) 

A 

where A is the domain of integration that is dependent both on the fluid 
inlet level and the cavitation boundary. 

Instantaneous load-carrying capacity is expressed as a ratio: 

(W)q > 0 represents the dimensionless steady state value of w evaluated 
for the same operating conditions and the same minimum film thickness. 

RESULTS AND DISCUSSION 

Results of numerical calculations that show the steady state performance 
of typical nonconformal contacts are presented in Tables 1 and 2, showing, 
respectively, the effects of the parameters Hq and a. 

The dynamic performance of hydrodynamically lubricated nonconformal con- 
tacts in combined rolling and normal motion is governed by the dimensionless 
normal velocity parameter q. It incorporates three major parameters, viz 
normal velocity, entraining velocity, and central film thickness. Representa- 
tive results for dynamic load ratio e, dynamic peak pressure ratio 5 , and 
the dimensionless locations of the peak pressure X and the film rupture 

r 

boundary are presented in Table 3 for a fixed value of dimensionless 
central film thickness Hq = l.OxlO"^ and with dimensionless normal veloc- 
ity parameter q varying between -1.0 and 0.75. In Tables 4 and 5 simi- 

-4 -5 

larly, Hq assumes the fixed values 1 x 10 and 1 x 10 , respectively. 

The dynamic load ratios of Table 4 are also displayed graphically in Fig. 3. 

The normal velocity parameter q clearly has a significant pressure-generating 
effect during normal approach and thereby increases the dynamic load ratio 3 
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with Increasing normal velocity parameter. On the contrary, 0 was signifi- 
cantly reduced during normal separation for q > 0. 

The magnitude of the peak pressure generated in a contact, its location, 
and the location of the film rupture boundary are also affected significantly 
by the normal velocity. Pressure distributions for various values of the nor- 
mal velocity parameter taken from Table 4 are shown in Fig. 4. The inlet 
meniscus boundary is not shown in these figures. During normal approach the 
film rupture boundary moves downstream into the exit region (or the divergent 
portion of the film) with reference to the minimum-film-thickness position. 
During separation it moves upstream, reaching the convergent portion of the 
film for higher values of normal velocity parameter. Similar observations 
were made by Dowson et al, (1976) fcr the lubrication of rigid cylinders in 
combined rolling and normal motion. The location of the pressure peak and the 
entire pressure distribution in the contact also shift accordingly. Thus 
superposition of normal motion on the entraining velocity alters both the mag- 
nitude and distribution of the pressure in the contact. 

Variation of peak pressure ratio c is shown in Fig. 5, again for data 
from Table 4. Pressures of the order of 3 to 4 times the corresponding peak 
pressure in the steady state situation are generated in the contact during 
normal approach. Relative reductions in the peak pressures of similar magni- 
tude occur during separation. 

The influence of the geometry parameter a on the dynamic load ratio 0 
and the peak pressure ratio ij is shown in Table 6 and in Figs. 6 and 7, 
respectively. Recall that for a > 1 the major axis of the film contours is 
transverse to the entrainment velocity, so for o » 1 the geometry resembles 
the usual case of a lubricated cylinder. During noriral approach e increases 
with increase in a up to a certain value and then iv tends to approach a 
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limiting value for further increase in a. Similarly, during separation, the 
e initially decreases with increasing a and then approaches a limiting value 
for higher values of a. Similar variations with x are observed for the 
dynamic peak pressure ratio c. Therefore it can be said that a signifi- 
cantly affects the dynamic performance of nonconformal contacts for values of 
a within the practically useful range 0.2 to 10. However, for a between 10 
and 35 the geometry parameter effects are small. For most practical situations 
the values of q and Hq fall within the ranges considered in the present 
calculations. For the sixty-three cases run e can be expressed as a function 
of q and a by a simple formula for the complete range of data presented in 
the Tables 3 to 6 as follows: 


e B (1.68 q) 


1/q 


for q 1 ^ 0 (11) 

Note that Eq. (11) does not explicitly contain the central film thickness 
parameter Hq for the following reason: While Hq appears both in the 
Reynolds equation and in the parameter q, comparison of Tables 3 to 5 reveals 
only a very minor explicit dependence of a upon Hq. For c this de- 
pendence is still less. It is therefore sufficiently accurate, except perhaps 
for larger negative q, to use values averaged with respect to Hq in de- 
riving a relationship such as presented in Eq. (11). 

The dynamic load ratio for rigid point contacts undergoing harmonic oscil- 
lations can also be determined by using Eq. (11). For small-ampl itude oscil- 
lations about the mean film thickness, the parameter q oscillates with 


amplitude q„,„ given by 

IhOX 


^max “ TIT 




uc 


( 12 ) 


mean 
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Since a negative q enhances the dynamic load ratio more than a positive q 
diminishes it, such oscillations vmt generate a higher average load-carrying 
capacity. 

It has been mentioned that significantly higher peak pressures can be 
generated in the contact during normal approach than during normal separation. 
Therefore the estimation of peak pressures in the contact is of practical im- 
portance for estimating the maximum stress in the contact and the fatigue life 
of the contact. Dynamic peak pressure ratio c is expressed by a simple 
formula as a function of the normal velocity parameter q and the geometry 
parameter a 

t s |a~0*032 (2q)|^^^ for q ^ 0.0 (13) 

Eqs. (11) and (13) are valid for the complete range of a between 0.2 and 35. 

Tables 7 to 0 show the percentage error between the values obtained from 
the preceding formulas and the computed values for both dynamic load ratio and 
peak pressure ratio. In general the agreement is excellent, and the formulas 
should find a good practical application. The slight explicit dependence of 
3 upon Hq has not been incorporated into the present formulas. 

The present formulas give excellent correlation for the dynamic peak 
pressure ratio c during normal approach. However, correlation does not 
appear good for higher values of q during separation. During separation for 
higher values of q (e.g., 0.5 and 0.75) the rupture boundary moves into the 
convergent portion of the fluid film and thereby shifts the pressure peak and 
the entire pressure distribution from the central film thickness into the inlet 
region. This alters the physical situation in the contact and perhaps warrants 
a closer look than was intended here. Rupture of the fluid film in the inlet 
region might lead to the failure of the lubrication process in the contact 
under certain situations. 


11 


f 


CONCLUSIONS 

A numerical solution to the hydrodynamic lubrication of nonconformal rigid 
contacts with an incompressible, isoviscous lubricant in combined rolling and 
normal motion has been obtained. The following conclusions were reached 
through a parametric study; 

1. Normal motion combined with pure rolling motion significantly in- 
creases load-carrying capacity peak pressure in the contact during normal 
approach. A correspondingly significant decrease in load-carrying capacity 
and peak pressure occurs during normal separation. 

2. The film rupture (or cavitation) boundary moves further into the exit 
region away from the minimum-film-thickness position during normal approach 
but moves toward the inlet and into the convergent portion of the film for 
higher normal vfJucities during separation. 

3. Incr 'fesing the geometry parameter increases dynamic load-carrying 
capacity and peak pressure during normal approach. Reverse effects are ob- 
served during separation. 

4. The dynamic load ratio a and the dynamic peak pressure ratio 5 
are very weakly dependent on central film thickness and can be adequately ex- 
pressed in terms of the dimensionless normal velocity parameter q and the 
geometry parameter o by simple equations as follows: 

a - ^ gg 1 /q 

t sech ( 2 q)M^^ 
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TABLE 1 . - STEADY STATE PERFORMANCE OF NONCONFORMAL CONTACTS - EFFECT 
OF MINIMUM FILM THICKNESS 


[Geometry parameter a ■ l.Oj inlet starvation parameter « 0.031; 
{liinensionless normal velocity parameter q 0.0 ] 


Dimensionless 
minimum film 
thickness , 

Ho 

Dimensionless 

load-carrying 

capacity, 

w 

Dimensionless 
peak pressure, 
Pmax 

Dimensionless 
location of 
peak pressure,^ 

Dimensionless 
location of film 
rupture boundary, ® 

Xr 

HASStSHil 

270.8 

3.77655xltf* 

-0.023 

+0.015 


1061.2 

1,2105x10° 

-.007 


1.0x10*3 

3493.76 

3 .8246x10' 

-.002 

+ .002 


3{-) indicates location tOo-ard Inlet side with respect to miniinuin-film- 
thickness position; (+) indicates location toward exit side. 


TABLE 2. - STEADY STATE PERFORMANCE OF NONCONFORMAL CONTACTS - 
EFFECT OF GEOMETRY 


[Dimensionless minimum film thickness Hn = 1.0x10"'*; inlet starvation 
parameter = 0.035; normal velocity parameter q a 0.0.] 


Geometry 

parameter, 

q 

Dimensionless 

load-carrying 

capacity, 

W 

Dimensionless 
peak pressure, 
Pmax 

Dimensionless 
location of 
peak pressure, 3 
Xp 

Dimensionless 
location of film 
rupture boundary,^ 
Xr 

0.2 

171.0 

0.4298x10® 

-0.008 

+0.003 

.4 

404.1 

.7173X10® 

-.008 


.6 

633.1 

.9217x10® 


)07 


.8 

S52.4 

1,0766x10® 



+.005 

1.0 

1061.2 

1.2105x10® 



+.005 

2.0 

1905.4 

1.6377x10® 



+.006 

5.0 

3660.0 

1.8479x10® 



+ .006 

10.0 

5646.6 

1,9958x10® 



+ .006 

20.0 

8168.3 

2.0498x10® 



+ .007 

35.0 

10987.2 

2.0820x10® 



+.007 


2{-) indicates location toward inlet side with respect to miniraum-film- 
thickness position; (+) indicates location toward exit side. 




















TABLE 3. - DYNAMIC PERFORMANCE OF NONCONFORMAL CONTACTS - EFFECT OF DIMENSIONLESS 

NORMAL VELOCITY PARAMETER 


[Dimensionless minimum film thickness Hn « 1,0x10"^; geometry parameter 
0 K 1»0; inlet starvation parameter Ti^f, 0,035.] 


Dimensionless 
normal velocity 
parameter,® 

Dynamic 
load ratio, 
e 

Dynamic peak 
pressure ratio, 
r, 

Dimensionless 
location of peak 
pressure,® 

^P 

Dimensionless 
location of film 
rupture boundary,® 
Xr 

-1.0 

2.98 

3.7408 

-0.008 


-.75 

2.316 

2.9258 

-.010 


-.5 

1.859 

2.2166 

-.013 


-.35 

1.403 

1.5531 

-.017 

+ .029 

-.1 

1.15 

1.2034 

-.020 

+ .021 

-.05 

1.073 

1.0984 

-.022 

+ .018 

+.05 

.931 

.9076 

-.025 

+.013 

+ .1 

.866 

.8223 

-.026 

+.010 

+.25 

.695 

.6161 

-.031 

+ .003 

+ .5 

.479 

.3460 

-.044 

-.010 

+.75 

.332 

.2017 

-.054 

-.020 


3(-} indicates location toward inlet side with respect to minimum-film- 
thickness position; (+) indicates location toward exit side. 


TABLE 4. - DYNAMIC PERFORMANCE OF NONCONFORMAL CONTACTS - EFFECT OF 
DIMENSIONLESS NORMAL VELOCITY PARAMETER 


[Dimensionless minimum film thickness Hn = 1.0x10"'’; geometry parameter 
a = 1.0; inlet starvation parameter = 0.035.] 


Dimensionless 
normal velocity 
parameter,® 
q 

Dynamic load 
ratio, 
e 

Dynamic peak 
pressure 
ratio, 

C 

Dimensionless 
location of peak 
pressure,® 

Xp 

Dimensionless 
location of 
film rupture 
boundary,® 

Xr 

-1.0 

2.783 

3.79 



-.75 

2.258 

3.0143 



-.5 

1.774 

2.2486 

O 

o 

I 

+ .014 

-.25 

1,348 

1.56 

-.005 

+ .009 

-.1 

1.129 

1.2053 


+ .008 

-.05 

1.063 

1.1002 

-.007 

+ .006 

+ .05 

.941 

.906 

-.008 

+ .004 

+ .1 

.885 

.8195 

-.008 

+ .003 

+ .25 

.738 

.5981 

-.010 

+.001 

+ .5 

.553 

.3462 

-.013 

-.003 

+ .75 

.422 

.2039 

-.017 

-.006 


®(-) indicates location toward inlot side with respect to miniinum-film- 
thickness position; {+) indicates location toward exit side. 
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TABLE B. - DYNAMIC PERFORMANCE OF NONCONFORMAL CONTACiS - EFFECT OF 
DIMENSIONLESS NORMAL VELOCITY PARAMETER 


tDimenslonless minimum film thickness Hq >• 1.0x10"^! geometry 
parameter a « X.O; inlet starvation parameter « 0.035,] 


Dimensionless 
normal velocity 
parameter, 
q 

Dynamic 

load 

ratio, 

B 

Dynamic 

peak 

pressure 

ratio, 

K 

Oimensionlesr 
location 
of peak 
pressure,® 
Xp 

Dimensionless 
location of 
film rupture 
boundary, ^ 
Xr 

-1,0 

2.62 

3.7588 

-Q.OOl 

+0,007 

-.75 

2.147 

2.9831 

-.001 

+ .006 

-.5 

1.71 

2.2184 

-.001 

+ .004 

-.25 

1.321 

1.5382 

-.002 

+ .003 

-.1 

1.12 

1.2111 

1 

+ .002 

-.05 

1,06 

1.1054 

1 

+ .002 

+.05 

.944 

.8984 

f 

+.001 

+ .1 

.894 

.8100 

-.003 

+.001 

+.25 

.754 

.6061 

-.003 

0 

+ .5 

.578 

.3540 

-.004 

-.001 

+ .75 

,454 

.2074 

-.005 

-.002 


indicates location toward inlet side with respect to 
min1mum-film-th1ckne,^s position; (+) indicates location toward 
exit side. 


TABLE 6. - DYNAMIC PERFORMANCE OF NONCONFORMAL CONTACTS - EFFECT OF GEOMETRY PARAMETER 


[Dimensionless minimum film thickness Hq = l.OxlO"'^; inlet starvation parameter Hin = 0.035.] 


Geometry 

parameter, 

a 




Dimensionless normal 

velocity parameter. 

a 




q - -1.0 


q 

= +0.75 





Dynamic 

peak 

pressure 

ratio, 

t 

Dimensionless 
location 
of peak 
pressure,® 
Xp 

Dimensionless 
location of 
film rupture 
boundary,® 

Xr 

Dynamic 

load 

ratio, 

B 

Dynamic 

peak 

pressure 

ratio, 

K 

Dimensionless 
location 
of peak 
pressure,® 
Xp 

Dimensionless 
location of 
film rupture 
boundary,® 

Xp 

0.2 

2.5971 

3.5697 

-0.003 

+0.019 

0.4485 

0.2232 

-0,017 


.4 

2.6684 

3.6646 



+ .021 

.4370 

.2148 




.6 

2.7195 




+ .023 

.4299 

,2101 



-0.006 

.8 

2.7533 

3.7704 



+ .024 

.4253 

.2063 





1.0 

2.7830 

3.7900 



+ .024 

.4220 

.2039 





2.0 

2.8594 

3.9227 

-.002 

+ .026 

.4121 

.1951 

-0.016 



5.0 

2.9380 




+.028 

.4032 

.1879 



-0.005 

10.0 

2.9770 

4.1242 



+ .029 

.3970 

.1840 





20.0 

2.9848 

4.1497 



+.029 

.3969 

.1829 





35,0 

2.9931 

4.1668 



+ .029 

.3966 

.1821 
























TABLE 7. - PERCENTAGE ERROR BETWEEN VALUES OF DYNAHIC LOAD RATIO AND PEAK 
PRESSURE RATIO GIVEN BY FORMULAS AND NUMERICALLY CALCULATED 
VALUES FOR VARIOUS DIMENSIONLESS NORMAL VELOCITY PARAMETERS 


Dimensionless 
normal 
velocity 
, parameter,® 

<1 


- 1.0 
-.75 
-.50 
-.25 
-.10 
-.05 
+ .05 
+.10 
+ .25 
+ .50 
+ .75 


Central film thickness, Hq 

1.0xl0~3 

1 . 0 x 10 -^ 

1.0x10-5 

1 . 0x10-3 

1 . 0 x 10 "'’ 

1.0x10-5 

Percentage of error in dynamic 

Percentage error in dynamic 


load ratio 

0 

pressure ratio c 

7.35 

-0.25 

-5.95 

0.57 

-0.73 


1.93 

4.55 

9.96 

6.93 

3.79 

4.88 

1.55 

6.43 

10.41 

7.42 

5.89 

7.3,3 

.43 

4.53 

6.66 

4.10 

3.64 

S.ll 

.07 

1.93 

2.75 

1.36 

1.20 

.72 

0 

.94 

1.23 

.60 

.44 

-.04 

.10 

-.96 

-1.28 

-.29 

-.11 

.73 

.34 

-1.81 

-2.80 

-.30 

.04 

1.21 

2.12 

-3.83 

-5.87 

.39 

3.43 

2.05 

-8.36 

-4.22 

10.58 

21.38 

21.31 

18.64 

21.70 

.47 

-7.1 

58.65 

56.76 

54.29 


H 


-) indicates location toward inlet side with respect to minitnum- 
film-thickness position; (+) indicates location toward exit side. 


TABLE 8 . - PERCENTAGE ERROR BETWEEN VALUES OF DYNAMIC LOAD RATIO AND 
PEAK PRESSURE RATIO GIVEN BY FORMULAS AND NUMERICALLY CALCULATED 


VALUES FOR VARIOUS DIMENSIONLESS GEOMETRY PARAMETERS 


Geometry 

parameter, 

Normal approach 
(q B - 1 . 0 ) 

Separation 
(q B 0.75) 

Normal approach 
(q = - 1 . 0 ) 

Separation 
(q = 0.75) 


Percentage error 
ratio 

in dynamic 

Percentage error in dynamic 
peak pressure ratio 5 

0.2 

2.18 

0.18 

0.10 

53.38 

.4 

1.41 

.31 

-.30 

54.74 

.6 

.61 

.41 

-.74 

55.49 

.8 

.21 

.60 

-.93 

56.42 

1.0 

-.25 

.38 

-.73 

56.76 

2.0 

- 1 .00 

.19 

-1.94 

59.06 

5.0 

-1.16 

- 1.02 

- 2.21 

58.82 

10.0 

-.54 

-2.09 

-1.80 

57.46 

20.0 

1.13 

-4.59 

-.22 

53.79 

35.0 

2.46 

-6.57 

1.17 

50.82 
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(a-2)x - 0 plane. 

(a) Two rigid solids separated by a lubricant film. 



(b-l)y-o plane (b-2) x - 0 plane. 


(b) Equivalent syst($m of a rigid solid near a plane separated by a lubricant film. 
Fig. Z - Contact geometry. 
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Fig. 3. - Variation of dynamic load ratio with dimension- 
less normal velocity parameter. Dimensionless central 
film thickness Hq-1.0x 10"^; dimensionless geometry 
parameter a- 1.0; and dimensionless inlet starvation 
parameter Hj^ ■ 0. 035. 
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Fig. 4. - Pressure distribution in 
thickness Hq- 1 . 0 x 10 "^; dimen 
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DIMENSIONLESS NORMAL VELOCITY PARAMETER, iql 

Fig. 5. > Variation of dynamic peak pressure ratio 
dimensionless normal velocity parameter. Dimension- 
less central film thickness FU - 1.0x10"^; dimension- 
less geometry parameter a - 1 . 0; and dimensionless 
inlet starvation parameter Hj^ - a 035. 
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Fig. 6, - Variation of dynamic load ratio with dimension- 
less geometry parameter Dimensionless central film 
thickness 1.0x10“^; and dimensionlecs inlet 
starvation parameter “ 0* 035. 



Fig. 7. - Var iation of dynamic peak pressure ratio with 
dimensioniess geometry parameter. Dimensionless 
central film thickness Hq “ 1.0x10"^; and dimen- 
sionless inlet starvation parameter Hj^ ■ 0.035. 
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